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RESUMEN 



ABSTRACT 

Expansion velocities (t4xp) of different ions and line widths at the base of 
the lines are measured and analyzed for 24 PNe with [WC]-type nuclei (WRPNe), 
9 PNe ionized by WELS (WLPNe) and 14 ordinary PNe. A comparative study 
of the kinematical behavior of the sample clearly demonstrates that WRPNe have 
in average 40-45% larger T4xp, and possibly more turbulence than WLPNe and 
ordinary PNe. WLPNe have velocity fields very much alike the ones of ordinary 
PNe, rather than the ones of WRPNe. All the samples (WRPNe, WLPNe and or- 
dinary PNe) show expansion velocities increasing with age indicators, for example 
<K;xp> is larger for low-density nebulae and also it is larger for nebulae around 
high-temperature stars. This age effect is much stronger for evolved WRPNe, sug- 
gesting that the [WC] winds have been accelerating the nebulae for a long time, 
while for non- WRPNe the acceleration seems to stop at some point when the star 
reaches a temperature of about 90000-100000 K. Non-WR nebulae reach a max- 
imum V^xp < 30 kms^^ while evolved WRPNe reach maximum 14xp of about 40 
kms~^. For all kind of objects (WRPNe and non- WRPNe) it is found that in aver- 
age I4xp(N"'") is slightly larger than V^xp(0+"'"), indicating that the nebulae present 
acceleration of the external shells. 

Key Words: ISM: KINEMATICS AND DYNAMICS PLANETARY 
NEBULAE: GENERAL PLANETARY NEBULAE: IN- 
DIVIDUAL: NGC6369, Ml-32, BD+30°3639, K2-16 — 
STARS: WOLF-RAYET 



1. INTRODUCTION 

Among galactic planetary nebulae with known 
central stars, less than 15% have nuclei presenting 
Wolf-Rayet features. All these stars have been clas- 
sified as belonging to the [WC] spectral type, show- 
ing almost pure helium and carbon in their atmo- 
spheres (e.g., Tylenda et al. 1993; Hamann 1997). 
The only possible [WN] central star in the Galaxy 
was reported by Morgan et al. (2003). 

Many recent studies have been devoted to these 
planetary nebulae (hereinafter WRPNe) and their 
central stars with interesting new results. Gorny & 
Stasihska (1995) found that the physical properties 
and chemical abundances of WRPNe were not dif- 
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ferent from those of non- WRPNe, except for higher 
expansion velocities. Since then, Peha et al.(1998, 
Paper I) found that central stars of different initial 
masses can pass through the same [WC] stage and 
Pena, Stasihska, & Medina (2001, Paper II) reported 
an unexpected behaviour of the [O III]/[N II] electron 
temperature ratios for WRPNe which does not seem 
to be present in non- WRPNe. In addition, Gorny 
et al. (2001) found that in an infrared diagram {H- 
K) versus (J-H), WRPNe spread over a larger zone 
than non- WRPNe which is probably more reflecting 
differences in dust and stellar properties than differ- 
ences in nebular properties. Acker et al. (2002) have 
concluded that WRPNe show similar expansion ve- 
locities but larger turbulent velocities that normal 
PNe, and recently Gesicki et al. (2003) concluded 
that the WRPNe tend to show strong acceleration 
at the ionization front and strong turbulent motions. 
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For the central stars, several scenarios to pro- 
duce such H-deficient low-mass stars have been pro- 
posed (e.g., Blocker 2001; Herwig 2001). In addi- 
tion, Gorny & Tylcnda (2000) and Dc Marco (2002) 
have provided further arguments for the existence of 
the evolutionary sequence: [WC]-late — > [WC]-early 
stars, ending with the PG 1159 type stars, proposed 
by Acker et al. (1996) and Hamann (1997). How- 
ever, this proposition is still debated (Hamann et 
al. 2005). Parthasarathy et al. (1998) have claimed 
that "weak emission line stars" (WELS), as defined 
by Tylenda et al. (1993), are an intermediate stage 
between [WC] and PG 1159 stars. However, Peha 
et al. (2001) have argued that PNe around WELS are 
more similar to ordinary PNe than to WRPNe and 
consider doubtful the proposition of Parthasarathy 
et al. 

Since 1995 we have performed systematic obser- 
vations of WRPNe, to gather a homogeneous high- 
resolution spectroscopic data set for studying the na- 
ture and evolutionary status of WRPNe. Prom these 
data, in Paper I we presented detailed photoioniza- 
tion models for a sample of very high excitation neb- 
ulae ionized by [WC 2-3] stars, and in Paper II we re- 
ported line intensities, physical conditions and chem- 
ical composition for 34 PNe with nuclei of all [WC] 
spectral types, including a few WELS. 

In this work, our data arc employed to analyze 
the kincmatical behaviour of 24 WRPNo in order 
to study the effects of the [WC] winds on the neb- 
ular shells. We analyze some phenomena affecting 
the nebular profiles such as expansion velocities, tur- 
bulence and high velocity material. In particular 
we address the problem of whether WRPNe present 
higher expansion velocities and/or turbulence than 
non- WRPNe. For our analysis we have included a 
sample of 14 ordinary planetary nebulae (i.e., whose 
nuclei are neither WR nor WELS) which are useful 
as a control sample to confirm that any special trend 
found for WRPNe is not an instrumental effect. Nine 
PNe ionized by a WELS (hereinafter WLPNe) have 
been included also in the analysis. 

In § 2, the observations and data reduction are 
described. The line profiles are discussed in § 3. The 
expansion velocities derived for the sample are pre- 
sented and analyzed in § 4, where a brief discussion 
regarding turbulence can be found. In § 5, we present 
the line widths at the base of lines, searching for the 
highest velocity material. The kinematical analysis 
of selected nebular lines is presented in § 6. Our main 
conclusions are summarized in § 7. Finally, two Ap- 
pendices were included where we present some ex- 
amples of our data and discuss some objects with 



particularly interesting line profiles. 

TABLE 1 

OBSERVING LOG.'^ 



PN G 


T Tan Q 1 n Q m o 


1 inG Hci'l'AG 

{<d<d/ mm/ yy j 


013.7- 


10.6 


YC 2-32 


02/11/00 


037.7- 


34.5 


NGC 7009 


01/11/00 


084.9 - 


03.4 


NGC 7027 


02/11/00 


103.7 + 


00.4 


M 2-52 


02/11/00 


104.4- 


01.6 


M 2-53 


01-02/11/00 


118.0- 


08.6 


Vy 1-1 


02/11/00 


130.3- 


11.7 


M 1-1 


01-02/11/00 


133.1- 


08.6 


M 1-2 


01-02/11/00 


147.4 - 


02.3 


M 1-4 


03/11/00 


159.0- 


15.1 


IC 351 


11/12/98 


166.1-1- 


10.4 


IC 2149 


13/12/98 


194.2 -1- 


02.5 


J 900 


02/11/00 


196.6- 


10.9 


NGC 2022 


02/11/00 


215.2- 


24.2 


IC 418 


12-13/12/98 


221.3- 


12.3 


IC 2165 


03/11/00 


243.8 - 


37.1 


PRTM 1 


12-13-14/12/98 


294.1-1- 


43.6 


NGC 4361 


04/10/99 



''Observing set-up: 2048 x 2048 (14 x 14 //m^) pixel CCD 
detector; slit dimensions: 4" x 13", oriented E-W; wave- 
length range 3360-7360 A, spectral resolution ~ 0.2 A. 
Extraction window: 4" x 3", for all the objects. 



2. OBSERVATIONS AND DATA REDUCTION 

High-resolution spectroscopic data were system- 
atically gathered with the 2.1 m telescope at the Ob- 
servatorio Astronomico Nacional, San Pedro Martir, 
B. C, Mexico on 1995 July 29-31, 1996 June 14-17, 
1997 August 3-4, 1998, December 11 14, 1999 Oc- 
tober 4-7, and 2000 November 1-3. The REOSC 
Echelle Spectrograph plus the University College 
London camera (described by Levine & Chakrabarty 
1993) were always used. Two different CCD detec- 
tors were employed: for observations on 1995 and 
1996 we used a CCD Tektronix of 1024x 1024 (23Atm) 
pixels, providing an average spectral resolution of 
0.3 A (about 18kms~^ in velocity) with a 150/um 
slit width, and since 1997 we used a CCD Thomson 
of 2048 x 2048 (14/im) pixels, giving a resolution of 
0.2 A (about llkms"^) with the same slit width. 

The observations for most of the objects pre- 
sented here, were described in Paper II, where a de- 
tailed list of objects, observing dates, exposure times 
and instrumental set-up were given. In Table 1 we 
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present the log of observations for objects not re- 
ported in Paper II. As usual we have taken at least 
two consecutive observations for each object. Differ- 
ent exposure times (from 2 to 15 minutes) were used 
to obtain a S/N better than 3 for the weak lines with- 
out saturating the most intense ones. The different 
exposures for the same object were averaged to in- 
crease the S/N. The slit was always E-W oriented. 
In this work, the reported data correspond always 
to those obtained for the central zone. For all the 
objects we covered a wide spectral range (at least 
from 3700 to 6800 A), therefore lines for high and 
low ionized species were obtained. This is important 
for studying the kinematics on different zones of the 
nebula where such species reside. 

A Th Ar lamp was used for wavelength calibra- 
tion in all the spectral ranges and standard stars 
from the list by Hamuy et al. (1992) were observed 
each night for flux calibration. Data reduction was 
performed using the package IRAF** and includes 
standard bias-subtraction and flat-field correction 
for all spectra. 

In Table 2, we present some nebular characteris- 
tics (diameter, electron temperatures from the [O III] 
lines, densities from the [S II] lines and ionic abun- 
dance ratios) for all the objects for which the kine- 
matical analysis is performed. Most of the data was 
taken from Papers I and II. For the new objects (Ta- 
ble 1) the nebular characteristics were derived by us- 
ing the same procedures as reported in those papers. 
In addition, in columns 8 to 11 we present the stellar 
classification, stellar temperature and wind parame- 
ters, taken from the literature. All these data will be 
used in the interpretation of the kinematical data. 

3. LINE PROFILES 

[WC] stars present intense and wide emission 
lines of He, C, and other elements, whose wave- 
lengths occasionally coincide with the wavelengths of 
some nebular lines. For instance, this is the case of 
the nebular He II, He I, and H Balmer lines which co- 
incide with stellar He II and He I lines, and especially 
it is the case of [Nil] A6583 and [S II] A6731 which 
appear on the top of some C II multiplet 2 lines. High 
spectral resolution is necessary to deblend the neb- 
ular and stellar components and to study the line 
profiles. This is essential for compact nebulae around 
[WC]-late stars where the stellar winds show termi- 
nal velocities of a few 100 kms~^, and the stellar 
lines present a FWHM of a few A. Therefore a spec- 
tral resolution better than 5000 is necessary. This 

''IRAF is distributed by NO AO, which is operated by 
AURA, Inc. under contract with the NSF. 



is one of the advantages of our spectra (resolution 
of about 20,000). In Appendix A, we present stellar 
and nebular lines for some of our objects to show 
that our spectral resolution allows us to safely de- 
blend both components for most of the objects. Thus 
the nebular lines have been measured without stellar 
contamination. 

Also our high spectral resolution, allows us to 
obtain well resolved nebular profiles for all the sam- 
ple, except for a couple of objects showing FWHM of 
lines smaller than our instrumental resolution. Thus, 
we have an extensive data set including different kind 
of objects, obtained in a consistent way which allows 
us to perform a reliable comparative analysis of the 
kinematics. First we will analyze the adequacy of 
deriving expansion velocities from the line profiles. 

3.1. Line Profiles and Velocity Field 
A nebular line profile depends on several inter- 
nal properties such as the nebular morphology, the 
density, thermal and ionization structures, the ex- 
pansion and turbulence velocities, etc. In addition, 
the observed profile also depends on the spectral res- 
olution and the position and dimensions of the slit: 
for instance, when the nebula is extended and the 
slit just covers the central zone, the classical double- 
peak profile is obtained. Many authors have derived 
the nebular expansion velocity, T4xp, as half the sep- 
aration of the peaks. On the other hand, when the 
slit includes the entire nebula, or we arc dealing with 
a filled nebula (with a very small inner hole) a single 
profile is usually observed. In such a case, the half 
width at half maximum (HWHM) has been used to 
measure I4xp- In our sample, many nebulae do not 
present a well defined shell, or they have knots at dif- 
ferent velocities, showing then complex profiles, for 
which is difficult to derive I4xp in a straightforward 
manner. 

An ad-hoc detailed kinematical model, includ- 
ing the density distribution (3D-morphology), and 
the ionization and thermal structures of the neb- 
ula, would be required to deconvolve the effect of 
each contribution to the line profile. An approach 
to this procedure, using a spherically symmetric 
model, can be found in Gesicki et al. (1996), Gesicki 
et al. (2003), and other papers by these authors. 
More recently, tomography techniques or 3D spatio- 
kinematic models have been used to calculate the 
line profiles of some individual nebulae (Sabbadin 
et al.2004 and references therein; Monteiro et al. 
2005). The main weaknesses of such a procedure 
are that the results are strongly model dependent, 
the derived model is probably not unique, and it is 
very much time consuming. 
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TABLE 2 

NEBULAR AND STELLAR PHYSICAL PARAMETERS^ 
PN G Name Te A^e He++/He+ 0++/0+ fli [WC]'' T. Voo log M Ref. 







(10* K) 


(cm-3) 






(arcsec) 




(103 K) 


(kms-l) 


(Mgyr-l) 




001.5 — 06.7 


SwSt 1 


1.06 


30000 


0.0 


0.3 


5 


9 


40 


900 


—6.72 


dMOl 


002.4 + 05.8 


NGC 6369 


1.01 


3220 


0.1 


15.7 


38 


4 


150 


1200 


—6.15 


KH97a 


002.2 — 09.4 


Cn 1-5 


1.09 


4700 


0.0 


82.8 


7 


4 


< 57 






TS94 


003.1 + 02.9 


Hb 4 


0.96 


4600 


0.1 


33.9 


7.2 


.3-4 


86 






TS94 


004.9 + 04.9 


M 1-25 


0.80 


10400 


0.0 


43.6 


3.2 


6 


60 






L96 


006.8 + 04.1 


M 3-15 


0.87 


7443 


0.0 


38.8 


4.5 


5 


55 






ZK93 


011.9 + 04.2 


M 1-32 


0.99 


6720 


0.0 


12.1 


9 


4-5 


66 






TS94 


012.2 + 04.9 


PM 1-188 


0.94 


200 


0.0 


0.03 


4 


10 


35 


360 


—5.70 


LH98 


017.9 — 04.8 


M 3-30 


1.03 


525 


0.9 


27.2 


20 


2 


97 






TS94 


027.6 + 04.2 


M 2-43 


1.05 


10100 


0.0 


14.3 


1.5 


8 


65 


850 


-6.08 


LH98 


029.2 - 05.9 


NGC 6751 


1.12 


2730 


0.0 


5.2 


20.5 


4 


135 


1600 


—6.12 


KH97b 


048.7 + 01.9 


He 2-429 


0.77 


6810 


0.0 


3.6 


4.2 


4-5 










061.4 — 09.5 


NGC 6905 


1.21 


1530 


2.9 


58.6 


40 


2-3 


141 


1800 


—6.32 


KH97a 


064.7 + 05.0 


BD+30°3639 


1.01 


21200 


0.0 


0.0 


7.7 


9 


47 


700 


-4.87 


L96 


068.3 — 02.7 


He 2-459 


2.00: 


17400 


0.0 




1.3 


8 


77 


1000 


-5.01 


L97 


089.0 + 00.3 


NGC 7026 


0.92 


3250 


0.2 


14.7 


20 


3 


130 


3500 


-6.34 


KH97a 


096.3 + 02.3 


K 3-61 


0.89 


1690 


0.0 




6.1 


4-5 










120.0 + 09.8 


NGC 40 


1.00 


1980 


0.0 


0.1 


48 


8 


78 


1000 


-5.62 


L96 


130.2 + 01.3 


IC 1747 


1.12 


2380 


0.2 


78.6 


13 


4 


126 


1800 


-6.58 


KH97a 


144.5 + 06.5 


NGC 1501 


1.05 


1020 


0.7 




52 


4 


135 


1800 


—6.28 


KH97a 


146.7 + 07.6 


M 4-18 


0.86 


6350 




0.0 


3.7 


10 


31 


160 


—6.05 


ciMC99 


161.2 — 14.8 


IC 2003 


1.28 


3840 


1.38 


21.4 


8.6 


3 


88 






TS94 


243.3 — 01.0 


NGC 2452 


1.30 


1590 


0.7 


6.0 


19 


2 


141 


3000 


-6.20 


KH97a 


352.9 + 11.4 


K 2-16 


1.17 


504 




0.4 


13.5 


11 


30 


300 


—6.36 


L97b 


009.4 — 05.0 


NGC 6629 


0.84 


3470 


0.0 


13.5 


15.5 


wl 


< 52 






TS94 


010.8 — 01.8 


NGC 6578 


0.79 


2270 


0.0 


240 


8.5 


wl 


65 






TS94 


011.7 — 00.6 


NGC 6567 


0.95 


4360 


0.01 


28.3 


7.6 


wl 


61 


1950 




TS94 


096.4 + 29.9 


NGC 6543 


0.79 


5640 


0.0 


75.4 


19.5 


wl 


< 66 


1900 


—7.4 


TS94 


100.6 — 05.4 


IC 5217 


1.12 


8850 


0.1 


93.9 


6.6 


wl 


72 






TS94 


159.0 — 15.1 


IC 351 


1.31 


2500 


1.0 


24.7 


7 


wl 


85 






TS94 


194.2 + 02.5 


J 900 


1.13 


1110 


0.7 


9.2 


9 


wl 


123 






TS94 


221.3 — 12.3 


IC 2165 


1.41 


3984 


0.8 


19.0 


9 


wl 


153 






TS94 


356.2 — 04.4 


Cn 2-1 


0.97 


5320 


0.07 


74.2 


2.4 


wl 


84 






P-M91 


013.7 - 10.6 


YC 2-32 


0.88 


3337 


0.03 




15 


pn 


68 






TS94 


037.7 — 34.5 


NGC 7009 


1.08 


4371 


0.3 


151.8 


28.5 


pn 


85 


2750 


—8.55 


TS94 


084.9 — 03.4 


NGC 7027 


1.64 


30000 


1.1 


2.8 


14 


pn 


175 






TS94 


103.7 + 00.4 


M 2-52 


1.41 


879 


1.3 


127.1 


14 


pn 










104.4-01.6 


M 2-53 


1.11 


496 


0.1 


1.6 


14.8 


pn 


112 






TS94 


118.0 - 08.6 


Vy 1-1 


0.99 


2101 


0.0 


56.6 


5.2 


pn 


32 






TS94 


130.3 - 11.7 


M 1-1 


2.7: 


3000 


46.7 


7.2 


6 


pn 


87 






TS94 


133.1 - 08.6 


M 1-2 




4264 


0.02 




18 


pn 


51 






TS94 


147.4 - 02.3 


M 1-4 


1.03 


6276 


0.06 




4 


pn 


67 






TS94 


166.1 + 10.4 


IC 2149 


0.95 


10000 


0.0 


6.8 


8.5 


pn 


< 49 


1290 


-7.7 


TS94 


196.6 - 10.9 


NGC 2022 


1.51 


800 


8.7 


62.5 


19 


pn 


114 






TS94 


215.2 - 24.2 


IC 418 


0.87 


18000 


> 0.1 


0.44 


12 


pn 


53 


1050 


-8.2 


TS94 


243.8 - 37.1 


PRTM 1 


1.63 


1200 


7.8 


> 300 


23 


pn 


90 






P90 


294.1 + 43.6 


NGC 4361 


1.93 


800 


> 19 


> 180 


63 


pn 


95 






TS94 



'^Nebular data are mostly from Paper II. Nebular diameter is from Acker et al. 1992. T,, Vc» and 
log M values were taken from literature, and the references are: dMOl: De Marco et al. (2001); dMC99: 
Do Marco & Crowthcr (1999); KH97a: Kocstcrkc & Hamann (1997a); KH97b: Koesterkc & Hamann 
(1997b); L96: Leuenhagen et al. (1996); L97a: Leuenhagen (1997a), L97b: Leuenhagen (1997b); LH98: 
Leuenhagen & Hamann (1998); P90: Pena et al. (1990); P-M91: Preite-Martinez et al. (1991); TS94: 
Tylenda & Stasinska (1994); ZK93: Zhang & Kwok (1993). 

''Central star type is as following: a number indicates the [WC] spectral type, wl is for WELS and pn 
for normal stars. 
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It is beyond the scope of this work to develop 
such a model for each of our objects. Our main aim 
is to derive reliable V^xp for different ions in order 
to compare the behavior of WRPNe with those of 
non-WRPNe and to address, if possible, the effects 
due to turbulence. 

To relate the line profiles with the values of the 

real expansion velocities for our objects in a con- 
sistent manner, we have used the code NEBU-3D, 
presented by Morisset et al. (2005), to compute a 
spatio-kincmatic model for analyzing the line profile 
behavior in a simple and common case: an ellipsoidal 
nebula ionized by a star with an effective tempera- 
ture of 80,000 K and a luminosity of lO'^"'' erg s^^. 
The chemical abundances are solar. The nebula has 
an internal elliptical hole with the major axis twice 
the minor axis. The orientation of the polar axis 
makes an angle of 45° with the line of sight (thus 
the nebula is in between "polar-on" and "face on" 
orientations). The density is constant in the radial 
direction, but varying with the polar angle according 
to N(H)in„er = N(Ho) X {Rq / Rinner? , where N(Ho) 
and i?o are the hydrogen density and the inner ra- 
dius at the equator. The dependence of N(H) with 
the polar angle is included in the parameter Rmner 
which follows the ellipsoidal form. 

Fig. 1 presents the H/? surface flux distribution, 
projected on the sky. The nebular dimensions pro- 
jected on the sky are about 5" x 7" (a distance 
of 1.87 kpc is assumed). Two "observing slits" (a 
3" X 4" slit, centered, and a similar one, off-center) 
have been superposed. These slits are used to extract 
"observational" data to compare with our own. 

The expansion velocity of the modeled nebula 
grows with the distance from the star as V = 60 
R/Rmax km s~^ {Rmax = 4.2 10^^ Cm is the dis- 
tance from the central star to the pole). A turbulent 
component of 0.5 km s~^ has been added to the ex- 
pansion. 

We consider that this hypothetical nebula repre- 
sents adequately the majority of nebulae. We have 
verified that for other morphologies (like spherical 
or bipolar nebulae) the results do not differ signifi- 
cantly from the conclusions presented here. A very 
simple spherical case has been analyzed by Gesicki 
& Zijlstra (2000) providing very similar results. An 
extensive catalog of profiles for elliptical and bipolar 
objects seen in different orientations and with differ- 
ent conditions will be published elsewhere (Morisset 
et al., in preparation). 

The main characteristics of our modeled nebula 
are presented in Fig. 2. The upper panels show the 
density, temperature and velocity structures as func- 



tion of the distance from the central star (in units of 
Rmax) along the equator (solid lines) and along the 
polar axis (dotted lines). The second row of pan- 
els presents the surface brightnesses computed for 
H/3, Hell A4686, and [Nil] A6583, as function of the 
distance from the central star along the same direc- 
tions. The third row presents the line profiles of H/3, 
HeIlA4686, and [Nil] A6583 as function of the ve- 
locity, obtained by integrating through the two slits 
shown in Fig. 1. The solid lines present the profiles 
through the centered slit (this slit has the same size 
and position as our extraction window), while the 
dotted lines show the profile through the off-center 
slit. The lowest row of panels shows the profiles 
through a 30" x 30" slit which includes the whole 
nebula. Superposed to these profiles we have traced 
horizontal lines indicating the values of 14xp as mea- 
sured with different methods: the heavy solid line is 
the "real velocity" computed as the average veloc- 
ity of the ion weighted by the emissivity of the line, 
the thin solid line is the peak-to-peak velocity, the 
dotted line is the HWHM velocity and the dashed 
line is the velocity at 1/10 the maximum intensity of 
the profile {Vio, see § 5). 

As it is observed in Fig. 2, our model shows that 
when the nebula is well centered in the slit and it is 
fully resolved (third row of panels) it presents sym- 
metric double-peak profiles and that half the peak-to 
peak separation is a very good approach for the "real 
velocity", underestimating it by less than 10%. The 
HWHM of the profile overestimates the "real veloc- 
ity" by about 10%. The profiles obtained when the 
slit is off-center (dotted profiles) are always asym- 
metrical and it is very difficult to measure the expan- 
sion velocities from them, as the separation between 
peaks is smaller than the centered case. 

When the nebula is fully included in the slit, 
it presents single Gaussian profiles and HWHM of 
the profiles is a good measurement of \4xp within 
20%. As expected due to the adopted velocity 
law, [N II] A6583 presents a higher expansion veloc- 
ity than HeIlA4686 and H/3, and the latter one 
shows much wider lines due to the thermal broad- 
ening (which is included). 

Vio is always larger than the 'real velocity' by 
about 30-40%, but it is closer to the velocity of the 
most external (although faint) zones and it is useful 
to detect high-velocity low-emission components in 
the gas. 

3.2. Our Observed Profiles 
Most of the profiles obtained for our objects fall 
into five different categories: a) symmetric double- 
peak profiles found mainly in extended objects, b) 
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Fig. 1. The H/3 surface brightness of the modeled nebula 
is shown. Two slits of 3" x 4" centered and off-center 
are superimposed. They are used to produce "observed 
profiles" (see § 3.1). 




Fig. 2. Upper panels show the density, temperature and 
velocity distributions in the modeled nebula. The sec- 
ond row of panels shows the surface brightness along the 
equator (solid lines) and the polar axis (dotted lines) for 
several emission lines. The third row of panels shows 
the line profiles obtained through the centered slit (solid 
lines) and the off-center slit (dotted profile). The hor- 
izontal lines represent the "real velocity" (heavy solid 
line), the peak-to-peak (thin solid line), the HWHM (dot- 
ted line) and Vio (dashed line) velocities, for the centered 
slit. The lowest panels show the line profiles through a 
slit of 30" X 30", including the whole nebula. 

single Gaussian-type profiles usually found for com- 
pact symmetrical nebulae, the prototype being M 4- 
18 whose [NII]A6583 profile is shown in Fig. 3a, c) 



single asymmetrical profiles as seen in Fig. 3b, usu- 
ally produced by asymmetric nebulae (or produced 
by aspherical nebulae where the slit did not cross ex- 
actly the center) , d) single profiles with wide wings or 
'shoulders' as shown in Fig. 3c (this could be due to 
high velocity gas or highly turbulent gas; high veloc- 
ity components have been reported for some objects, 
e.g., De Marco et al. 1997), e) complex profiles show- 
ing several components which usually correspond to 
knotty nebulae as the one presented in Fig. 3d. 

For the first two categories, 14xp can be mea- 
sured safely in the classical form: from half the peak- 
to-peak separation for the wcU-rcsolved objects and 
from the HWHM for objects showing single lines 
with no asymmetries nor complex profiles. Our mea- 
surements arc presented in the next section. Objects 
with very asymmetrical profiles, for which it is very 
difficult to define Vexp in a simple way, were excluded 
from the T^xp analysis. For all the objects, the type 
of line profile is indicated in the last column of Ta- 
ble 5. 

For the analysis, the profiles of H/3, Hell A4686, 
[O III] A5007, and [N 11] A6583 lines were measured 
when available. The purpose for selecting these 
lines was to study possible kinematical differences 
between lines arising from ions located at different 
distances from the central star. For instance, while 
the H/3 line is mapping the whole nebula. He II A4686 
and [O III] A5007 lines give information about the in- 
ner and intermediate regions, and [N II] A6583, of the 
outer regions. 

4. EXPANSION VELOCITIES 
4.1. Objects with Split Profiles 

Twelve objects of our sample (six WRPNe, two 
WLPNe, and four ordinary PNe) present symmetri- 
cal split profiles. According to our model such pro- 
files are observed for well resolved nebulae and half 
the peak-to-peak separation is a good approach for 
Kxp- Values for [O III] A5007, H/?, [NII]A6583, and 
HeIlA4686 lines are presented in Table 3. Typical 
errors are ±5 kms"-'^. The [WC] type of the cen- 
tral stars or an indication if the star is a WELS or a 
"normal" central star is listed in Col. 3. 

The line widths (FWHM) for each component 
of H/3, AVi, and AV2 (in kms~-^), are also listed. 
The instrumental and thermal widths have been 
subtracted from each line width by assuming that 
they add in quadrature. The instrumental widths 
were measured from the comparison lamp lines and 
they resulted to be 18 kms"'^ for observations prior 
to 1997 (23/x pixel size CCD) and 12 kms~^ since 
1997 (14/i pixel size CCD). The thermal contribu- 
tion to line width, for any ion can be expressed as 
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AVth = ^lA^tc/A (kms where te is the electron 
temperature in units of 10^ K and A is the atomic 
mass (Lang 1980), therefore for H/?, AVth = 21. 4^/*^ 
while for [OIII], AVth = 5.35Vt^, etc. The electron 
temperature employed for each object has been listed 
in Table 2. 

The individual line widths in Table 3, in the 
range of 18-30 kms~^, are therefore mainly due 
to turbulence and possible velocity gradients within 
the shell. In this sense notice that the model pre- 
sented in § 3, which has a velocity gradient pro- 
ducing AV > 20 kms~^ from side to side of the 
shell, pre diets a FWHM of about 5 6 km s^^ for the 
double-peak profiles (third row of panels in Fig. 2), 
after subtracting the thermal width. Our objects in 
Table 3 present much higher AV, possibly produced 
by turbulence. 

Table 3 (see also § 6.1) shows that, for a given ob- 
ject, the expansion velocities from different ions are, 
in general, very similar, although most of the objects 
have slightly larger Vexp([NII]) than Vexp([OlII]) or 
14xp(H/3) which, as in our model, might be indicat- 
ing that the expansion velocity increases with the 
distance from the central star. This fact, predicted 
by hydrodynamical models as due to acceleration of 
the external nebular material, was already noticed 
by Wilson (1950) for his sample of PNe. 

Our Vexp values are equal, within uncertainties 
with those from the catalogue by Weinberger (1989), 
which have been included in Col. 10 of Table 3. In 
Col. 11 we have also included the expansion velocities 
found by Acker et al. (2002), based upon modeling 
the velocity field, for the few objects in common; 
their Vexp values and ours coincide, except for the 
faint and extended WRPN K2-16 for which Acker 
et al. give two different values. K2-16 is discussed 
in the Appendix B. 

For the sample in Table 3, we found that WRPNe 
have VcKp values (as measured from [O III] lines) 
ranging from 24 to 44 kms~^ with an average of 
36 kms^^, while WLPNe and ordinary PNe show 
Vixp from 17 to 26 kms^^ with an average of 21.5 
kms-i. That is < Voxp(WRPNe)> is 67% larger 
than <Vcxp(non-WRPNe>). Therefore, definitely, 
WRPNe in this sample have larger expansion veloci- 
ties than non- WRPNe (WLPNe and ordinary PNe). 
Also a systematic trend of higher V^xp with earlier 
[WC] type seems to be present. This result will be 
strengthened in the following sections. 

In addition to higher V^xp. WRPNe seem also 
to show larger line widths than non-WR objects. 
This is particularly true for WRPNe with early 
[WC] stars such as NGC 6905, NGC 6369, NGC 6751, 



and NGC 1501, which in principle are more evolved, 
probably indicating large turbulence in these nebu- 
lae. 

In conclusion, WRPNe with split profiles are 
showing higher Vexp and probably more turbulence 
than WLPNe and ordinary PNe. This does not 
agree with the results presented by Acker et al. 
(2002), who have concluded that the expansion ve- 
locities in WRPNe and non- WRPNe are similar 
but WRPNe arc much more turbulent than non- 
WRPNe. Also their models would indicate that or- 
dinary PNe present acceleration in the outer nebular 
zones, while WRPNe do not, which is not the case 
for our objects. 

4.2. Objects with Single Profiles 

The FWHM of [O III] A5007, H/3, [NII]A6583, 
and HeIlA4686 were measured for all the objects 
showing single symmetrical profiles ([N II] A6548 was 
measured when the stellar C II lines were severely 
blended with [NII]A6583. See Appendix A). Instru- 
mental and thermal widths were subtracted as indi- 
cated in § 4.1. We have assumed that HWHM rep- 
resents the expansion velocity plus turbulence. The 
results are presented in Table 4 where, in Col. 8, we 
have included V!xp([OlII]) from the compilation by 
Weinberger (1989) or from more recent literature and 
in Col. 9 the values reported by Acker et al.(2002). 

Although in this case it is more difficult to com- 
pare kinematical data measured by different authors 
using different procedures, it is found that the ma- 
jority of our values coincides within a factor of 1.2 
with those from the literature. For several objects 
(Ml-32, He 2-459, Kl-61, and YC2-32), Vexp are 
determined for the first time. 

Typical errors for our measurements are ±6-8 
kms~^ for [Nil] and Hell lines and ±4 kms~^ for 
[O III] and up lines. We have searched for system- 
atic errors (as, e.g., a dependence of Vsxp as a func- 
tion of electron temperatures) and found that these 
errors are negligible. Also possible artificial broaden- 
ing of nebular lines due to contamination with stellar 
lines, is discarded as a source of important system- 
atic errors (see Appendix A). Systematic trends of 
V^xp with electron density and stellar temperature 
are found for all the sample. These are real trends 
and they are discussed in the next sections. 

Vexp averages of the 14 WRPNe in Table 4 are 23, 
23, and 25kms-\ for [Olll], H/3 and [Nil], respec- 
tively. The corresponding averages for the 5 WLPNe 
are 19, 18, and 22 and for the ordinary PNe are: 19, 
17, and 20. WLPNe and ordinary PNe show a very 
similar behavior and can be put together in a single 
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Fig. 3. Position-velocity diagrams and integrated profiles from the extracted spectra are shown featuring the different 

types of line profiles found in our sample. The different typos arc: a) simple Gaussian [Nil] A6583 profile of M4-18, 
b) asymmetrical [N II] A6583 profile of Hb 4, c) [O III] A5007 high velocity wings found in M 1-32, d) knotty structure 
and complex [O III] A5007 profile in NGC 2452. At the top of each figure, a velocity scale helps to visualize the velocity 
profile. The corresponds to the maximum line intensity. 



sample, non-WRPNe, showing averages of 19, 18, 
and 20kms-i for [OIII], H/3, and [Nil]. Therefore, 
the averages for WRPNe are 20-25% larger than the 
corresponding averages for non-WRPNc indicating 
that WRPNe have wider lines. In this case we can- 
not disentangle expansion velocities from turbulence, 
but this result is consistent with the results found 
in § 4.1 of higher expansion velocities and probably 
higher turbulence in WRPNe than in non-WRPNe. 

As in the case of double-peak objects, here the 
V^xp average for [Nil] is about 2kms~^ larger than 
the corresponding values for [O III] , regardless if they 
are WRPNe or non-WRPNe, confirming that the ex- 
ternal shells might be expanding faster in all kind of 
nebulae. 



5. LINE WIDTHS AT THE BASE OF LINES 

Considering that an important fraction of PNc 
shows kinematical features like ansae, BRETS, 
FLIERS, jets, high velocity ejecta, etc., (see Bal- 
ick & Adam 2002 for a definition of these terms), all 
of them produced by non-spherical, bipolar, multi- 
polar or asymmetrical ejections of the central star, 
it seems likely to expect some additional perturba- 
tions in the velocity field of WRPNe as compared 
with non-WRPNe, caused by the [WC] wind. For in- 
stance, the hydrodynamical models of Garcia-Segura 
&: McLow (1995) show that the expansion of a hot 
bubble pushed by a WR wind results in a filamentary 
broken shell with large turbulence. 

Therefore, to perform a kinematical analysis of 
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TABLE 3 

EXPANSION VELOCITIES AND TURBULENCE FOR NEBULAE WITH DOUBLE PEAK LINES 



PN G 




Name 


[WC]'' 


Vcxp (k 


ms ^) 


b 






14,p(0III)^ 


cxp 










fo ml 


H/3 


[Nil] 


He II 


H/? 


up 


W89 
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HAi A no 
UDi.4 — uy 





inVjU oyuo 
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A A 
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iNljO DoDy 


A 
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36 
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O / 


QQ 

oo 


Z4 


QQ 
00 
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4Z 




noo o r\K 
uzy.z — Uo 


o 
y 


iNljiVj D/Oi 


A 
4 


42 


40 


/in 

4U 




z4 


zy 


/in 

4U 


/I 1 


144.0 + UD 


c 





4 


40 


40 


A 
4z 


OO 


on 


ly 


/ 


/in 


120.0 + 09 


8 


NGC40 


8 


30 


26 


26 




18 


< 15 


29 




352.9+ 11 


4 


K2-16 


11 


24 


24 


26 




22 


16 




34 or 38 


100.6 - 05 


4 


IC5217 


wl 


18 


18 


23 


18 


16 


20 


23 




011.7-00 


6 


NGC6567 


wl 


17 


18 


27 




< 20 


< 20 


19 




037.7- 34 


5 


NGC 7009 


pn 


20 


18 


21 


16 


< 15 


< 15 


21 




196.6 - 10 


9 


NGC 2022 


pn 


24 


21 


21 


23 


16 


20 


27 




243.8 - 37 


1 


PRTMl 


pn 


24 


23 




23 


28 


26 






294.1 + 4.S 


(i 


XGC4.S()1 


])U 


2() 


22 




22 


■■V2 


20 


21 





''Central star type as in Table 2. 

''Uncertainties: ±5kms~'^ for Kxp, and ±8kms~^ for line widths. AV^s were corrected for 
thermal and instrumental widths. 

'^References are: W89: Weinberger (1989); A02: Acker et al. (2002). 



all our sample, including all these effects, we decided 
to consider the nebular gas showing the highest ve- 
locity relative to the star. This can only be measured 
at the base of the line profiles but to avoid any noise 
disturbance due to the low signal at the base of the 
lines, we measured such a velocity as half the width 
of the line profile at one tenth maximum intensity 
(Vio), after subtracting the instrumental and ther- 
mal widths as described in § 3.1. A similar treat- 
ment was previously used by Dopita et al. (1985) to 
determine expansion velocities of a sample of PNe in 
the SMC. 

Vio was determined for all the objects, including 
those cases where a single Gaussian was adequately 
reproducing the line profiles or when the lines showed 
a double peak. In the latter case, half the full width 
at 1/10 intensity of each component was added to the 
peak-to-peak separation. It can be easily calculated 
that, for a single Gaussian profile, half the full width 
at 1/10 I max can be expressed as: 

Fio = 0.911 FWHM (kms-^); 

while for double peak profiles (two Gaussians), it is 
found that: 

Vio = Kxp + 0.455 (AFi + AV2) (kms-i) . 

We have used these expressions to derive Vio for our 
objects in Tables 3 and 4. For non-Gaussian, wing- 



extended, or complex profiles, Vio was measured di- 
rectly from the line profiles. 

Fio derived in such a way includes not only Vexp 
of the shell, but also the turbulence of the gas and 
faint high velocity components if they exist. The 
values are listed in Table 5. Uncertainties are about 
±5 kms~^ and larger if marked with a colon. In the 
last column of Tabic 5 we describe the line profiles: 
s for single, d for double peak, c for complex, h for 
high velocity components. 

6. KINEMATICAL ANALYSIS 

Many authors have searched for correlations 
among the nebular velocity field and the nebular 
and stellar parameters. The reported results are 
not always consistent. For instance, Sabbadin (1984) 
and others reported a correlation between the neb- 
ular radius and V^xp, which has been interpreted as 
changes in Vc^p as a consequence of nebiilar cvohi- 
tion. Gesicki & Zijlstra (2000) have not found such 
a correlation for a sample of bulge PNe with well de- 
termined distances. Also acceleration of the exter- 
nal shells (appearing as larger V^xp for N+ and 0+, 
than for 0++) seems to be a well established fact 
for the majority of PNe. Nevertheless, from their 
kinematical model. Acker et al. (2002) reported no 
evidence of such an acceleration for WRPNe, but re- 
cently Gesicki et al. (2003) — using a similar model 
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TABLE 4 

EXPANSION VELOCITIES IN OBJECTS WITH SINGLE SYMMETRICAL LINES 
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NGC 7027 


pn 


22 


21 


28 


22 


22 




103.7 + 00.4 


M 2-52 


pn 


19 


19 


16 


22 


8 




104.4-01.6 


M 2-53 


pn 


17 


21 


20: 


22: 


11 




118.0-08.6 


Vy 1-1 


pn 


12 


13 


18: 




10 




130.3-11.7 


Ml-1 


pn 


31 


30 


22: 




39 




147.4-02.3 


M 1-4 


pn 


15 


14 




11 


14 




166.1 + 10.4 


IC 2149 


pn 


< 9 


11 


18 




< 10 




215.2 - 24.2 


IC 418 


pn 


< 9 


8 


14 




< 6 


15 



''Central star type as in Table 2. 

'^Typical errors arc +6 8 krns^^ for [Nil] and Ho II, and ±4 kms""' for H/3, and 
[OIII]. Errors are larger for the cases marked with a colon. 

"^Data are from Weinberger (1989) unless indicated otherwise. Other references are: 
1: Dc Marco et al. (2001); 2: De Marco & Crowther (1999); 3: Gesicki & Zijlstra 
(2000). 

'^A02: Acker et al. (2002). 

" Brycc & Mellema (1999) reported V'exp([0 III])=36, 14xp([N II])=28 kms"^ for 
BD+30°3639. 



for a sample of 14 PNe fomid that the line pro- 
files of most nebulae in their sample can be better 
reproduced by adopting an U-shape velocity field, 
whether the central star is a [WC] or not. 

Here we use our consistent data to study the be- 



haviour of Vcxp vs. several properties of the nebu- 
lae and their central stars, searching for the physical 
causes affecting the velocity field in WRPNe. 
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Fig. 4. Kxp from H/3 are compared to T4xp from: [Nil], 
He II, and [O III]. WRPNe are marked with black circles, 
WLPNe are marked with open circles and ordinary PNe 
arc marked with open triangles. Least squares fits are 
shown as solid lines. Typical error bars are shown in each 
case. The WRPN conspicuously out of the Vixp([0 III]) 
versus Voxp(H5) relation is BD+30°3639, discussed in 
Appendix B. 

6.1. Vcxp from Different Ions 

In Fig. 4, the values of Ve^p derived from [N II] , 
He II, and [O III] lines are plotted against the value 
derived from H/3. Different kind of objects (WRPNe, 



WLPNc, and PNc) are shown with difFcront symbols 
and typical error bars are presented in each graph 
(±8 kms"^ for [Nil] and Hell, and ±4 kms"^ for 
[O III] and H/3) . An immediate result from this fig- 
ure is that Vexp for WRPNe extends over a wide 
range from lOkms"^ up to 50 kms~^, while l^xp 
for ordinary PNc and WLPNc arc always lower than 
30 kms~^. For the whole sample, the averages of 
Vexp for WRPNe are about 40-45% larger than the 
averages found for non- WRPNe. Notice that, as al- 
ready said, (Kxp (WRPNe)) / ( Kxp (non- WRPNe)) 
is larger for double-peak objects than for objects 
with single profiles. As double peak nebulae are 
more extended and presumably older objects, this 
would indicate that evolved WRPNe expand faster 
that evolved non- WRPNe. 

Fig. 4 shows that Kxp for [Nil], Hell, and 
O III follow a very close 45° linear correlation with 
Kxp (H/3), although in the case of [O III] versus H/3 
there is one WRPN, BD-F30°3639, that conspicu- 
ously departs from this behavior (see more details 
in the Appendix B). Interestingly, we do not find 
significant differences among WRPNe, WLPNe, and 
PNe; all kind of objects present the same correlations 
within uncertainties. The least square fits shown in 
Fig. 4 are: 

• Kxp([NlI]) = (0.92 ±0.08) Kxp(H/3) + (4.0 ± 
2.0), 

r = 0.94 , 

. Kxp(Hell) = (0.89 ±0.12) Kxp (H/3) + (3.0± 
2.9), 

r = 0.95 , 

• Kxp([OlII]) = (1.02 ± 0.08) Kxp(H/3) + (0.0± 
2.0), 

r = 0.98 . 

In Fig. 4 and from the least square fits it can 
be seen that the [Nil] vs. H/3 graph shows a 
larger dispersion but, within the errors, we find 
that in average, Kxp ([Nil]) is slightly larger than 
Kxp (H/3). [Hell] vs. H/3 correlation also might 
indicate that (Kxp (He II)) is slightly larger than 
(Kxp (H/3)) but the uncertainties arc larger in this 
graph and the result is not conclusive. On the other 
hand, Kxp([0lll]) is equal, within uncertainties, to 
Kxp (H/3). As already said, higher Kxp for low ion- 
ization species were found in other PN samples (Wil- 
son 1950; Gesicki et al. 2003), and could be due to ac- 
celeration of the external zones relative to the inner 
shells or to higher turbulence of the external zones 
due to shocks or expansion of the ionization front 
against the neutral external material. We found that 
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such a phenomenon is occurring in WRPNe as well 
as in non- WRPNe. 

The same statistical results are obtained if we use 
Vio as given in Tabic 5. Vio averages for WRPNe are 
also about 40-45% larger than the respective aver- 
ages of non- WRPNe. As Vio is computed for a larger 
number of objects the results are strengthened. Nev- 
ertheless the main importance in the analysis of Vio 
compared to Kxp is the possible detection of very ex- 
tended wings or high velocity components produced 
by low emission gas, which are only noticeable when 
Vio is employed. This is discussed in more detail in 
the next sections. 



6.2. Correlation o/V^xp with the Nebular Density 

In Fig. 5 (top), electron densities derived from 
the [S II] line ratio are plotted vs. V^xp from H/3. The 
graph is largely scattered. A wide range of velocities 
corresponds to a given density, but it is evident that, 
regardless of the stellar type (WR or not) the nebular 
density and V^xp are anti-correlated. As high den- 
sity PNe are considered young nebulae, the behavior 
found in Fig. 5 (top) is confirming an evolutionary ef- 
fect in the sense that low density (therefore evolved) 
nebulae expand faster and possibly have larger tur- 
bulence than high density (young) objects. 

In this figure, young high-density WRPNe (all 
of them ionized by [WC]-late stars), present a wide 
range in velocities (from 12 to 30 kms~^) while most 
young ordinary PNe and WLPNe are constrained to 
V^xp < 20 kms"'^. Certainly, the powerful [WC] wind 
seems to largely affect the velocity fields in seemingly 
young PNe. For the objects lying in the low density 
zone (log Uc < 3.5), WRPNe appear much faster 
than normal PNe. Our sample includes only three 
WLPNe in this zone, but they are showing the same 
behavior as ordinary PNe. 

Fig. 5 (bottom) shows the behavior of density 
against Vio (H/3). The velocity ranges are wider in 
this case, with Vio from 20 to 65kms~-'^ for WRPNe, 
while all WLPNe and PNe, except two, show values 
of Vio < 45kms~^. In both figures, the upper open 
triangle at log = 4.48, corresponds to the young 
PN NGC 7027 which shows an extraordinarily large 



V^xp for its density. It could be indicating the exis- 
tence of a dense stellar wind. This was also proposed 
by Keyes & AUer (1990) in order to obtain an en- 
hanced stellar radiation field shortward 130 A, nec- 
essary to fit a photoionization model for the nebula. 
In recent observations of the central star, Hubeny et 
al. (1994) did not detect such a wind, however its 
presence in the near past cannot be discarded. 
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Fig. 5. Electron densities derived from [S II] line ratio 
are plotted versus Vexp (top) and Vio (bottom) from H/J. 
Symbols are as described in Fig. 4. Despite the large 
scatter, both quantities are anticorrelated. The upper 
open triangle at log ne=4.48 in both panels, corresponds 
to the young PN NGC 7027. 

6.3. Correlations o/V^xp with Characteristics of 
the Central Stars 

6.3.1. Stellar Temperature and Nebular Expansion 

The fundamental stellar parameters of many 
[WC] stars have been calculated from non-LTE ex- 
panding atmosphere models by several authors (see 
references in Table 2). It is found that stellar tem- 
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peraturcs, T*, arc closely related to [WC] types 
(e.g., Koesterke 2001). It should be noticed that 

given by models corresponds to the temperature 
at the base of the dense wind (optical depth r ^ 20, 
e.g., Koesterke & Hamann 1997b) and these are not 
necessarily the temperatures of the ionizing fluxes 
"seen" by the nebulae. In fact there arc some well 
documented [WC] stars where the temperature at 
r = 2/3 is much lower than and, as a result, 
the PNe show very low excitation. A typical case is 
NGC6751 which does not present nebular He^^ al- 
though it is ionized by a [WC 4] star of T*=135,000K 
at T = 20. Its T, at T = 2/3 is, however, as low 
as 55,500 K (Koesterke & Hamann 1997b). This is 
very different from what occurs in ordinary optically- 
thick non-WRPNc, where the He II Zanstra tempera- 
ture and the stellar effective temperature are similar 
(Tylenda & Stasinska 1994). 

To study the behaviour of Vc^p and Vio func- 
tion of T*, we gathered the values of T* from the 
literature (see Table 2). For [WC] stars which have 
been modeled, we adopted the model values. For the 
non-modeled [WC] stars and the other objects we 
have listed the He II Zanstra temperatures from the 
same data base as compiited by Tylcmda & Stasinska 
(1994). Although this is a good approximation for 
the eff'ective temperature for non-WRPNe, it could 
be only a lower limit for the stellar temperature of 
[WC] stars. 

Fig. 6 presents the behaviour of V^xp versus 
(top) and Vio versus (bottom). In spite of the un- 
certainties affecting , a clear tendency is noticed in 
both graphs showing that hotter stars ionize nebulae 
with larger velocities and there is no object with a 
large T* showing V^xp < 20 kms~^. The scatter is 
however very large and the correlation coeflacient for 
these graphs are 0.64 for 14xp versus T* and 0.55 for 
Vio versus T,. 

All nebulae, regardless if they are WRPNc or not, 
display this behaviour. Two main groups can be dis- 
tinguished in the graphs: In the first, with below 
100,000 K, some WRPNe show velocities as low as 
non-WRPNe and both type of objects appear mixed 
in the lowest velocity zone; in the the second group, 
with > 100,000 K, WRPNe and non-WRPNe are 
well separated and WRPNe fall on a zone of much 
higher velocities than non-WRPNe. 

Assuming that the stellar temperature is an indi- 
cator of the age of the central star (which is the case 
for central stars evolving at approximately constant 
luminosity), it is clear that the older the star, the 
faster the nebular expansion, indicating an acceler- 
ation of the nebula with age for all type of stars. 



However for non-WR objects the curve in Fig. 6 
seems to flatten when the stars reach ~ 100,000 
K. That is, non-WRPNe would stop accelerating 
and their Vcxp reach an upper limit of about 25-30 
kms~^ while the acceleration continues for WRPNe. 
Evolved WRPNe show I^xp of about 40 kms'^ 

This apparently continuous acceleration from 
young to evolved WRPNe is what one might ex- 
pect if [WC] stars evolve from late to early types. 
A detailed hydrodynamical evolutionary model, in- 
cluding the kinematical effects of a long-term [WC] 
wind, should be compared with the behavior de- 
tected in Fig. 6, allowing us to confirm or discard 
the proposed evolutionary trend. Such a model is in 
progress (Medina et al., in preparation) and will be 
published elsewhere. 

The larger dispersion found in the Vio versus T» 
graph shows that some objects, in particular some 
young WRPNe and also a few ordinary PNe with 
relatively cold stars, present large turbulence and/or 
high velocity components since early phases of evo- 
lution. These phenomena are present in many PNe 
irrespective of their age. 

Acker et al. (2002) have not found any correlation 
between I^xp and in their data, because they have 
used the stellar temperatures obtained from pho- 
toionization models. That is, they used "seen" 
by the nebulae, which, as explained above, could be 
several thousands degrees lower than the tempera- 
ture at the base of the wind of the [WC] star. Any 
possible relation is thus destroyed. 

6.3.2. [WC] Mass-Loss Rate and Wind Velocity 
versus Nebular Characteristics 

It is worth to notice that the kinetic energy car- 
ried by the particles in the stellar winds can easily 
ionize the innermost gaseous material thus perturb- 
ing the ionization equilibrium. Therefore, to ana- 
lyze this possibility, in Fig. 7 we present the neb- 
ular ionization degree He++/He+ of WRPNe as a 
function of the mass-loss rate, M, and the termi- 
nal velocity of the wind, Vqo (these parameters are 
listed in Table 2). It is evident in these graphs that 
highly ionized WRPNe are associated with stars with 
low mass- loss rates (log M < — 6) and Vqo larger 
than 1800 km s~^. This is expected if we consider 
that such wind parameters are typical of [WC]-early 
(therefore hot) stars, which produce very excited 
nebulae with large I4xp- On the other hand, [WCj- 
latc stars which have higher M and lower Vqo only 
produce low ionization PNe with no traces of He"'""'". 
Therefore, we have found that the powerful stellar 
winds do not contribute significantly to the ioniza- 
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Fig. 6. Vcxp (top) and Vio (bottom) are plotted versus the 
stellar temperatures. Symbols are as described in Fig. 4. 
Despite the large scatter, velocities and temperatures are 
correlated. The correlations coefficients are 0.64 for the 
top figure and 0.55 for the bottom figure. 

tion of nebulae and the ionization can be safely at- 
tributed to the UV stcUar photons. 



Fig. 7. Nebular ionization degree (He^'^'/He"'') is plotted 
against some stellar parameters such as: a) the mass- 
loss rate (M), and b) the wind terminal velocity (voo)- 
Symbols are as described in Fig. 4. 



6.4. The Mechanical Energy of the Stellar Wind 
and the Expansion Velocities 

In Fig. 8 we plot Vcxp vs. the rate of mechanical 
energy of the stellar wind, L^, calculated as 1/2 M 
(see data in Table 2). A very clear tendency is 
found showing that stars with large Lyj have nebulae 
with large Kxp- The linear fit has a correlation co- 
efficient R = 0.71, indicating that, despite the scat- 
ter, the trend is real and much of the acceleration 
of the nebulae around [WC] stars should be due to 
the wind mechanical energy. Unfortunately only one 
WLPN and three ordinary PNe of our sample have 
published M and Voo, nevertheless, as expected con- 



sidering their low 
zone. 



all of them lie in the low 



7. MAIN CONCLUSIONS 

In this work we have analyzed high spectral res- 
olution data to study the kinematical behavior of a 
large sample of planetary nebulae around [WC] stars, 
as well as PNe around WELS and ordinary PNe. The 
main conclusions are as follows: 

• From a spatio-kinematic model representing an 
ellipsoidal nebula we have established that Vc^p 



34.5 35 35.5 

Log L_w (erg/s) 



Fig. 8. Vixp is plotted against the log of the rate of me- 
chanical energy of the stellar wind, L^. Symbols are the 
same as in Fig. 4. It is evident the relation between both 
parameters. The linear fit is X4xp=7.8(log L^) —249.3 
kms~^, with a correlation coefficient R = 0.71, therefore 
the large \4xp of WRPNe can be attributed to accelera- 
tion by the [WC] stellar wind. 



can be measured as half the peak-to-peak sepa- 
ration of nebular lines for well resolved nebulae 
observed through a centered slit. The HWHM 
for objects with single symmetrical profiles is 
also a good indicator of Vcxp, in particular for 
the cases of low turbulence. 
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Considering the above, we measiircd T4xp of sev- 
eral ions, for the objects in our sample showing 
double-peak or single symmetrical profiles. For 
a number of objects (K2-16, PRTMl, Ml-32, 
M2-43, He 2-429, He 2-459, K3-61, NGC6578, 
and YC 2-32) these values are presented for the 
first time. 

For all kind of objects (WRPNe and non- 

WRPNc) it is found that the velocities of the 
different ions are very similar but in general; 

Vexp([Nll]) > yexp([0lll]), Fexp(H^). 

This indicates that all kind of nebulae presents 
acceleration of the external shells. 

We have shown that WRPNe have a different 
kinematical behavior than non- WRPNe. For 
WRPNe Vexp extends over a wide range from 
10 kms~^ up to almost 50 kms~^, while for or- 
dinary PNe and WLPNe Vcxp is always lower 
than 30 kms~^. Average expansion velocities, 
measured from the peak-to-peak separation of 
split lines or from the FWHM of lines, are 40- 
45% larger in WRPNe than in non- WRPNe. 
This difference is larger for objects showing 
double-peak profiles. In addition to larger Kxp, 
WRPNe seem to present more turbulence than 
ordinary PNe. These results are based mainly 
on our sample with double-peak profiles for 
which T4xp and turbulence effects are easily sep- 
arated. 

It is worth to remark that PNe around WELS 

were found to have a velocity field very much 
alike the one of ordinary PNe, rather than 
the one of WRPNe. This evidence shows 
that WELS, contrary to the suggestion by 
Parthasarathy et al. (1998), do not belong to 
the [WC] family. 

Searching for evidence of high velocity material 
and perturbations of the velocity field, half the 
widths at 1/10 of maximum intensity of the neb- 
ular lines (Vio) were measured. We have found 
that (Vio) for WRPNe are also about 40-45% 
larger than the averages for non- WRPNe. But 
in particular Vio is large for some young com- 
pact WRPNe and a few ordinary PNe which 
would have large turbulence or high velocity ma- 
terial in the central zones. This phenomenon is 
present in many planetary nebulae irrespective 
of their evolutionary age. 



• Strong evidence indicating that the nebular ex- 
pansion increases with the nebular age (in- 
dicated by the nebular electron density), re- 
gardless of the nebular type (WRPNe or non- 
WRPNe) has been found. 

• The same is found when analyzing the tem- 
perature of stars: the older the central star 
(as measured by its temperature), the faster 
the nebular expansion. This behavior occurs 
in non- WRPNe as well, but these objects seem 
to accelerate to a maximum Vexp of about 25- 
30 kms~^ when their stars reach T* of about 
90,000-100,000 K, while WRPNe continue ac- 
celerating up to reach T4xp~ 40 kms~^. The 
acceleration of nebulae is closely related to the 
mechanical energy of the stellar winds. [WC] 
winds accelerate more and apparently for longer 
time their nebTilac. This constitutes an indirect 
evidence supporting the validity of the proposed 
evolutionary path: [WC]-late stars [WC]- 
early stars. Another evidence for this was es- 
tablished years ago through the determination 
of the chemical abundances in [WC] stellar at- 
mospheres. The results showed that C abun- 
dances were similar in [WC]-early and [WC]-late 
stars (De Marco 2002). However those results 
have been recently questioned by Hamann et 
al. (2005) on the base of improved non-LTE 
atmosphere models. These models show lower 
C abundance in apparently evolved [WC] -early 
stars than in less evolved [WC]-late ones, while 
one should expect the opposite. A further revi- 
sion of this subject is certainly worthy. 
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APPENDIX A 

A. Deblending Nebular from Stellar Lines 

In Fig. 9 we present observed profiles of Ha, 
and [Nil] AA6548,6583 lines, and [Sll] AA6717,6731 
for different [WC]-type objects. Some stellar lines 
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(usually very wide) appear at the same wavelength 
that the nebular lines and could affect the nebular 
line measurements if the spectral resolution is not 
adequate. These figures show that our resolution 
of about 18,000, allows to safely deblend the stel- 
lar and nebular lines, even in the cases of [WC]-late 
stars whose stellar winds show velocities of a few 100 
kms~^ and thus their lines have FWHM of a few A. 
In a few extreme cases, where [N II] A6583 appears 
severely blended with the stellar C II lines, as the 
case of M4-18 or NGC40, the line width of [Nil] 
A6548 was measured. 

APPENDIX B 
B. PROFILE OF INDIVIDUAL OBJECTS 

B.l. The Peculiar M 1-32 and BD+3(P3639 

The WRPNe M 1-32 and BD-F30°3639 show sin- 
gle symmetrical profiles, but they present V\q from 
[O III] larger than Vio from H/3 by about SOkms"-'^. 
For BD-F30°3639 this is also the case for Kxp- It 
is well known that this object shows this anomaly 
(Bryce & McUema 1999 and references therein) and 
we have found that M 1-32 presents a similar behav- 
ior. 

The spatially resolved study of BD-h30°3639 by 
Bryce & MeUema (1999) shows that the N+ shell is 
more spatially extended but less spectrally extended 
than the O"'""'" shell. The [N II] velocity ellipses pre- 
sented by these authors appear to be almost open- 
ended, with [O III] emission emerging from the gaps. 
Bryce & Mellema measured the expansion veloci- 
ties from [Nil] and [OIII] profiles finding 28±1 and 
36±1 kms~^, respectively, which coincide within un- 
certainties with our values of 23 and 43kms~"'^ for 
the same ions. Acker et al. (2002) reported 14xp of 
27kms~^ with a turbulent velocity of 15kms~^ for 
this object. 

M 1-32 has very unusual profiles characterized 
by a very-narrow intense component (with FWHM 
smaller than 25kms~'^) and very wide wings, with a 
full width at the base of about 125kms~^ (Fig- 3c). 
The wings are much wider in [O III] A5007 and 
He I A5876 than in H/? and they should be produced 
by high velocity gas in the unresolved central zone. 
Turbulence also can produce large wings, but a huge 
turbulence velocity should be required to create the 
wide wings observed (Morisset, in preparation). 

Both objects are low ionization nebulae 
(BD+30°3639 in particular), where 0++ lies 
in the inner nebular zone which is probably interact- 
ing with (and being disturbed by) the [WC] wind. 
It seems evident that the mechanical energy of the 



wind is strongly affecting the inner zone in both 
nebulae. 

B.2. Objects Showing Single Asymm,etrical Profiles 

Fig. 3b shows the single asymmetrical profile 
of [N II] A6583 for Hb 4. This is a high excitation 
Type I planetary nebula with a bright nebular core 
and outer low ionization ansae (Corradi et al. 1996; 
Gongalves et al. 2001). Lopez et al. (1997) showed 
that the ansae are produced by coUimated outflows 
with radial velocities of ±150 kms^^, relative to the 
central core. They measured an expansion velocity 
of 21.5 kms^^ for the core. The model by Acker 
et al. (2002) indicates expansion and turbulence ve- 
locities of 16 and 14 kms^^, respectively. Our l^xp 
value for H/3 is 16 kms~^. 

Our position- velocity diagram for the central core 
(Fig. 3b) shows a structure similar to a broken or 
incomplete shell, with a bright compact knot in the 
red side and fainter extensions to the blue. This 
structure is more evident at low-ionization species 
and produces asymmetrical single line profiles with 
a Vio of about 32 kms~^. 

Another nebula with similar asymmetrical pro- 
files is IC 1747, where we also detected a bright knot 
and faint extensions to the blue. Extracted profiles 
show an intense component with an extended blue 
wing. Images by Balick (1987) show a knotty ring 
nebula. 

The ordinary PNe Ml-2, M2-53, and IC2149 
show similar single asymmetrical profiles. 

B.3. Objects Showing High Velocity Extensions 
or Very High Turbulence 

The most conspicuous object presenting high ve- 
locity extensions is M 1-32, described above. Other 
objects showing the same phenomenon are: He2- 
459, although the profiles are not as wide as in M 1- 
32, and M3-15, where two very faint high- velocity 
knots (located in the central zone) at —90 and 
80 kms~^ are detected. In the latter case, the knots 
are faint with intensities below -j^/max! thus, only the 
velocity of the bright component is reported in Ta- 
bles 4 and 5. M 3-15 is also peculiar in the sense that 
this nebula shows an extraordinarily low Ne abun- 
dance (Paper II). 

There are three of ordinary PNe showing very 
large Fiq. They are PRTMl, NGC4361, and Ml- 
1. The first two are considered halo PNe, and are 
both density bounded nebulae. The low ionization 
species are not detected in these objects (Pena et al. 
1990; Torres-Peimbert et al. 1990). Vazquez et al. 
(1999) reported high velocity material and a com- 
plex kinematics in NGC 4361. Similarly M 1-1 shows 
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Fig. 9. Profiles of Ha, and [Nil] AA6548,6583 lines, and [S II] AA6717.6731, showing the nebular and stellar lines 
for different [WC]-type objects ([WC4] NGC6369, [WC8] NGC40, [WC 10] M4-18. The nebular lines can be safely 
deblended from the wide stellar lines even in the cases of slow stellar winds as NGC40 and M4-18, except for [Nil] 
A6583. [N II] A6548 can be measured instead. 



a very high excitation with faint emission of the low 
ionized species. It is usual to find fast T4xp in den- 
sity bounded objects (Acker et al. 2002). In these 
thin low-mass nebulae the ionization front has bro- 
ken through the outer layers which apparently ex- 
pand fast having no neutral material around. 

B.5. K2-16: a [WC]-Late Object 
with Double-Peak Profiles 

This is a very-faint low-ionization PN with a shell 
morphology of 13" in diameter. The shell structure 
as well as the low electron density are unusual for a 
nebula around a [WC 11] star. In this sense K 2-16 is 
similar to the also extended and low-density PM 1- 
188, ionized by a [WCIO] star (Pefia 2005). K2-16 
is the only nebula around a [WC]-latc star showing 
double peak profiles. The star is very bright (V = 
12.75 mag from Acker et al. 1992) and it presents 
intense emission lines severely contaminating some 
nebular lines, specially in the zone of [Nil] A6583 
and [Sll] A6731. Therefore we took special care in 
measuring the nebular line widths in a zone at 4-5" 
outside the stellar emission. In Fig. 10 we present the 
extracted spectra showing "the nebular only" and 
"the nebular + stellar" emissions. 



Our V^xp value of 24 kms"-"^ is much smaller than 
the expansion velocity determined by Acker et al. 
(2002), who proposed two models for this object: 
one with V^xp = 34 kms~^ and a turbulent ve- 
locity of 12 kms~^ or alternatively, a model with 
Vexp = 38 kms~^ and an acceleration producing a 
difference of 35 kms""*^ between the velocities of the 
inner and outer radii. It seems that the stellar lines 
affecting the nebular [N II] A6583 has introduced an 
additional uncertainty in Acker et al. model. 

The kinematics of K 2-16 could be explained by 
an old thin shell strongly accelerated by the [WC] 
wind of this apparently low-mass slowly-evolving 
star. A very simple calculation to support this idea 
can be made by assuming that the nebular mass of 
K 2-16 is about 0.1 Mq. Then the mechanical energy 
in the shell would be about 10"^^ '^^ erg. Thus the 
wind (with mechanical luminosity 1/^=10^^ '^^ erg 
s~-^ as calculated from Table 2) should have been 
blowing for about 1500 yr in order to provide the 
energy in the shell (an efficiency of a 100% is as- 
sumed). A similar calculation for M4-18, gives an 
age of about 580 yr. In conclusion K2-16 should be 
several times older than M4-18, for possessing such 
a low density, large diameter and high I4xp- 
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Fig. 10. Nebular only (up) and nebular + stellar (low) ex- 
tracted spectra in the zone of Ha for K 2-16. The nebular 
spectrum was extracted at about 4-5" from the central 
star. The nebular lines show double-peak. The stellar 
-I- nebular spectrum shows the intense stellar lines C II 
lines well mixed with the nebular emission, preventing to 
measure [Nil] A6583. 

B.6. Peculiar Line Profiles in NGC 6369 
This extended nebula shows a pecuhar position- 
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velocity structure (see Fig. 11). 
well resolved picture of this object 



beautiful 
found at 

http://heritage.stsci.edu/2002/25/index.html The 
nebula consists of a faint central zone, a very bright 
extended ring, faint filamentary extensions and ansa- 
type structures that were studied by Gongalves et al. 
(2001). Monteiro et al. (2004) modeled this object as 
a clumpy hourglass shaped nebula. In the HST pic- 
ture, the ionization structure is easily seen showing 
that the low ionization species are located mainly in 
the outer zone. 

In Fig. 11 we show that, at the central zone 
(0", where the central star is located) well split lines 
for all the ions are detected, with V^xp of about 
36kms^^. At about 3-4" to the East from the cen- 
tral star we detect a very bright zone (corresponding 
to the inner side of the ring) emitting mainly in high 
excitation lines (HI, [OIII], He I, and Hell). These 
lines present single profiles. The low ionization lines 
( [N II] , [O II] , [S II] are weak here and they present 
double peak. This should be the zone where Acker 
et al. (2002) obtained their spectrum with Ha show- 
ing a single component and [N II] A6583, a complex 
profile (see their Fig. 2). Due to the complexity of 
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Fig. 11. Position- velocity contours and the extracted 
spectra from (a) [O III] A5007 and (b) [N II] A6583 are 
shown. The position of the central star is at 0" in both 
cases. In (a), the solid line corresponds to the emission of 
the bright knot at 3" East from the central star, whereas 
the dashed double-peak line corresponds to the emission 
at the central star position (0"). This emission is shown 
5x its actual intensity to feature both lines in a single 
plot. In (b), both emissions at 0" and 3" show the same 
split-line spectrum. The bright knot is not detected in 
[N II] A6583. 

the profile in this zone they could not model this 
object. 
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TABLE 5 
Vio FOR THE SELECTED IONS 
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''Central star type as in Table 2. 

''Uncertainties of Fio values are in average about of 4-5 kms"'^, and larger for the cases marked 
with a colon. 

'^Line profiles are as following: s: single compact, a: asymmetrical, h: high velocity extensions, 
d: double peak, and c: complex. 



